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ABSTRACT 


The ability to control a viscoelastic field is an essential concept that defines some 
important restrictions and potentials of the influenced material. This thesis investigates 
the controllability of three popular constitutive models under homogeneous extensional 
and shear flows via the Lie bracket method. The constitutive models are as follows: the 
Phan-Thien-Tanner model; the Johnson-Segalman model; and the Doi model. The effect 
of extensional flow on these models and the effect of shear flow on the Doi model have 
not been explored previous to this work. The main contribution of this thesis is to 
characterize the submanifolds in the state space on which the non-Newtonian flow fields 
are weakly controllable. This kind of approach based on the control Lie algebra can be 


applied to a wider variety of complex models. 
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I. INTRODUCTION 


A. BACKGROUND AND PURPOSE 


Controllability of viscoelastic fluids is a characteristic that is significant in the 
design of desired materials. Though the theory of control for linear systems has been 
aptly developed, nonlinear problems remain unfamiliar. Recent studies from Renardy [1] 
investigated the nonlinear case and characterized controllability and the associated 
reachable set of several popular constitutive models under homogeneous shear flow of 
viscoelastic fluids. Analysis revealed that for some of the models the states in stress space 
are specified by a positive definiteness inequality of the stress tensor and are reachable 
from a given initial condition. Renardy followed up with additional research which 
extended to the control of nonhomogeneous shear flow, explicitly to the upper convected 


Maxwell fluid [2]. 


The purpose of this thesis is to broaden Renardy’s research [1] on the 
controllability of viscoelastic fluids by incorporating a variety of homogeneous steady 
flow fields. This investigation encompasses nonlinear characteristics, thus the method of 
analysis presented is principally based upon the nonlinear geometric control theory, 
which differs from Renardy’s work [1]. An additional difference is the definition of 
controllability, where this paper focuses on a local adaptation of the concept, known as 
weak controllability. The definition of weak controllability has been extensively used in 
nonlinear control theory and engineering applications. In practice engineers prefer to 
attain a distance target by making a sequence of local movements to reduce the side 
effects of model uncertainties, system perturbations and sensor noises. So, weak 


controllability provides more practical information for engineering applications. 


Though the study here can be applied to other complicated constitutive models, 
this thesis will focus on the following three representative fluids: the Phan-Thien- 
Tanner; the Johnson-Segalman; and the Doi. The first two are popular models for 
viscoelastic fluids whereas the last one is a classical model for liquid crystalline 


polymers. 


B. CONTROLLABILITY OF CONTROL SYSTEMS 


For reader’s convenience, we present some basic definitions and theorems 


adopted from [4]. More details can be found in [3]. 


Let 
x= f(x)+>)9,(x)u, (1) 


be a general nonlinear control system, described to be affine in control where x is the 














state variable, and the control variables are u,¢U for i=1,...,.m. Overall, the state 
variable obtains a value on a n-dimensional smooth manifold denoted by M. For 
systems described by (1), controllability is a significant property, and plays a crucial role 
in many control problems. Controllability distinguishes the manipulation capability of 
the system under control, and the stabilization of unstable systems via feedback. There 
exists a substantial amount of literature on controllability theory. For this thesis, the 
geometric approach from Isidori [3] will be adopted. Definitions are as follows: 
Definition 1 [4]. A point x, in M is characterized as reachable from x, if there 
exists piecewise constant input functions, u, =a,(t), such that the corresponding 
trajectory x(t) with initial state x(0) =x, reaches x(T) =x, in finite time for some 


O<T. 


For nonlinear control systems, the global reachability is usually very difficult to 


prove. Instead, a practical solution is to study weak controllability. 


Definition 2 [4]. A control system is characterized as weakly controllable within 


some open subset SCM, if for each x, in S, there is an open region U,, of x, such 
that the set of points reachable from x, to x(t) within U,, encloses at least one open 


subset of M. 


A weakly controllable system indicates that the locally reachable states create a 


“solid”, non empty region. The theory of controllability is more restrictive than weak 


2 


controllability such that it mandates any two points in M must be reachable from each 
other. In linear control systems, concepts of controllability and weak controllability are 
interchangeable. Alternatively, with nonlinear control systems (such as the ones 
presented in this thesis), characterizing the controllability of the system would require 
information on the global geometric properties of the vector fields. For general control 
systems, this remains to be a challenging quandary. Interestingly enough, by exploiting 
the dimension of the control Lie algebra and a distribution produced by the vector fields 
related to the control system, the weak controllability can be characterized. This brings 


to light several significant properties of a control system [4]. 

Consider the notion of smooth vector fields defined on some manifold M . These 
vectors, such as f(x)and g,(x) of (1), can be given different algebraic structures which 
can be useful in understanding the controllability of the system. Lie Algebra, for instance 


is such a structure in which the product of vector fields is defined by their Lie bracket, 


[f.g]. Itis within this operation where the minimum subalgebra containing /, g,....,2,, 


is described as the Control Lie Algebra, symbolized by @. Thus, within manifold M , at 
every point x, the vectors in @ span a vector space. These vectors produce a distribution 


defined by A, (x)= span{X (x)|X is a vector field in 2 }. 


Definition 3 [4]. Let n be the dimension of manifold M. A control system 


satisfies the controllability rank condition (CRC) on an open set SCM _ if 


dim(A, (x))=n for allx in S. 


Theorem 1 [3]. Let n be the dimension of manifold M . A sufficient condition 
for control system (1) to be weakly controllable on an open set S is for the 


dim(A 3 (x)) =n for allx in S. Or said differently, the control system (1) is weakly 


controllable if it satisfies the controllability rank condition on S. 


C. FLOW FIELDS 


Figure | and Figure 2 are graphs of the flow fields that will be used in this thesis. 


Figure 1 represents a homogeneous extensional flow where the elongational rate is 


3 


constant and the velocity is described as: v -(104-100 2), Upon investigating 


Figure 1, one can see that there is a vertical contraction and a horizontal extension. 
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Figure 1. An Example of An Extensional Flow 


Figure 2 represents a homogeneous shear flow where shear rate is constant and 


the velocity is described as: v = (7(t) y,0). In investigating Figure 2, one can see that 
there is a shearing motion in the horizontal direction. 


Shear flow and extensional flow have many physical applications. For example, 
in extrusion manufacturing the flow away from boundaries is well approximated locally 
by elongation; in film and sheet manufacturing, or in squeezing flow between parallel 
disks, there are interior flow regions approximated by elongation. In fact, many 
manufacturing process flows can be decomposed as a combination of shear flow and 


extensional flow. 
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An Example of A Shear Flow 


-0.5 





THIS PAGE INTENTIONALLY LEFT BLANK 


I. THE PHAN-THIEN-TANNER MODEL 


Of the nonlinear viscoelastic constitutive equations, one of the most applied 
differential type is the Phan-Thien-Tanner model [5]. The Phan-Thien-Tanner model has 
two parameters that control non-linearity and has the ability to fit, to some extent, the 
shear and elongational properties independently. It gives, however, spurious oscillations 


during start-up of shear flow. Thus the model is still liable for improvement. 


The Phan-Thien-Tanner model is described by the following state equation: 
T -(Vv)T -T (Vv) + AT +x(trT )T =2uD (2) 


where the definition of its element are as follows: 


T. 


Lie 
T : Stress tensor > T -| - ‘| (where 7,, =T7,, due to symmetry) 
21 22 


T : Derivative of the stress tensor with respect to time 


v: Velocity where v =(¥,,v,) 


Bie Si 
. ‘ Ox Oy 
Vv: Gradient of the velocity > Vv = 
Die HO 
Ox oy 


(Vv 4 : Gradient of the velocity transposed 
A: Relaxation rate 

«x: Transposed velocity gradient tensor 
trT : Trace of the stress tensor 

i: Elastic modulus 


Vv +(Vv y’ 


D: Rate of deformation tensor (the symmetric part of Vv ) > D= 


y: Control input (closely related to velocity) 


A. HOMOGENOUS EXTENSIONAL FLOW 


A fluid in a homogeneous extensional flow with rate y(t) is defined by velocity 


y(t)} 1 O 
v= [n3-702) , where the velocity gradient tensor is Vy = nae al . Recall, 


a schematic plot of extensional flow is given in Figure 1 of section I. This velocity is 


applied to the system (2), whereby the system’s components are investigated. 


nr. 0() Spl, *)eoservrmely 2 


HO Ti, | A 4) E ' k ‘] 
=>T +——+ + +(A+K«(trT = uy |(t 
2 | Ti. T,, qT, TT, ( ( ) T, T,, ( ) O: 2=1 
-7(t) 21; 0 T, Tp . 1 0 
+ A T,,+T. a 
2 0 21, Heese) Ty Hy 4 
To S02] 2.0) 
=>T =-(A+K(T,,+T. - a t | 4 ey 
aust Nr Ts 0-1] [0 -Z,| 


T, T Ee oF M+T, 0 |. 
=| n]—(asx(n eta) ac rr ro GB) 


12 22 








The system’s components are as follows: 
=f =—-(A+ «(T+ T))T +(H tT, )7(t) 
=>T, =-(A+«(T, +T))T -(44+Tx)7(t) 
=>, =-(A+«(f,+T,))T 


For mathematical convenience, the following notion is introduced: 


T,, —(A+«(x,+x,))x, M+X, 

= -(A+«(x,+x,))x, a —(#+%,) y(t) 
—(A+«(x,+%,)}%, 0 

The following Lie Brackets are computed: Fag and |F[ 4-8 ]]- 

First Lie bracket: | fg )=Ve-f-V -g 
i & 3G —(A+«(x,4+%,))x, 

=| #58 |= O -1 O}- —(A+K«(x, +x) X, 
0 0 


0 —(A+K«(x, +2) 


—(A+«(2x,+x,)) -KX, 0 U+tx, 
= -KX, —(1+«(x, +2x,)) 0 ‘| —s- X, 


—KX, —-KX, —(A+K(x, +x,)) 


—_— 
ba 
ay 


: —(A+ K(x, +%)) x —(A+«(2x, + x,))(4+™)—«x% (-w-%) 
=| 8 |= (A+ (x, +%,))x, - —Kx, (ut+x,)-(A+«(x, +2x,))(-4-%) 
0 -Kx,(+x,)—Kx, (-u-%,) 


-(A+K(x, +Xx,))x,+(A+« (2x, +X,) (+x, )+ Kx, (-u-x,) 
=| fig | = (A+ K(x, +%,))x, + Kx, (ut x,)+(A+«(x, +2x,))(-w- 2) 


KX, (+x, )+KX,(—u—- x) 
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Kx, (% —x,)+ (A+ Kx, a) 
(4) 


o[ral-fstioa}ensenos 


K(x, —X,)X, 


Second Lie bracket: lf f8]| =V| f.8|-f-Vi | 7.8 


=|f[f@ 


—Kx, — (A +Kx, + Kx,) + «xx, |-|-(A+K(x, + X_)) % 


—(A+K(x, +X,)) x; 


Kx, + (At KX, + KX, )— KX, —(A+«(x,+x,))x, 
ley 
K(%,—%) x, 


—(A+«(x,+x,))%, [ Kx, + M(A+ Kx, + KX) —KXX 
-V| -(A+x«(x,+%,)) x |-|-«x) -u(AtKx, +x, )+«xx, 


—(A+ K(x, +%,)) x, K(x, — xy), 


UK + KX, —2KX, — UK + KX, 0 





cca LK — KX, 0 | [-(4+« he 


KX, -KX; Kits) 


KX, (A+ x(x, +2x,)) 0 KX, (%,—x,)— (A+ Kx, + Kx) 


es KX, 0 f e oahetionces) 


ea KX, (A+(x, +x) KX3(% — Xp) 


o[7{tall-stento M+ X,)x, +( - an X, 


Hy et (Hay 


(2x, +uU-Xx,)x,+(u-%) . | 
Kx, (% — 2) [ 2 +2 ( x, +.x,) ]+ u( A+K(x,+%)) 
+| «x, (x, —x,) )[ A+ 2«(x, +x) ]- U(A+K(x,+ x))) 


KX; (x —X, )fA + 2K (x, + %)'| 


2 (4) ela) 
eke lf z|| =-K(A+K(x, vat = = em ) 


2%;( %=3;,) 


KX, (xX —x,)[A+2K(x, +x,)|+u(A+«(x, +x,)) 


( 
+] (4-9) [44 24( +4) ]-a(4 «(4 45)) 
KX, (x, —x,)[A+2«(x, +x,)| 
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X, (4, — 3) + 4(% +%) 
=[7[Pa]}--eaentverf lana 
Kx, (4, +.)(4- 4) + H(AtK(% oa), 
+] 1X (% +35 )(%4 = %,)—M(A+ K(x, +x))) 
Kua +5) ip) 


nie(A44(5;445))[ (4-4) a} nn (+2)(55—¥) 


> EAERal = -K(A+ K(x, ra) (x,-x,)+ a}en’s (x, + x, )(%,—%) 


—K (A+ K(x, +X,)) x, (% —x,) +x, (% +x,)(% —%) 


Als ALE 


| -KAx, (x, -x,)+wA(A+K(x, +x,)) 
=| 7 f.2]|- —KAx, (x,-x,)- MA(A +(x, +%)) (5) 
| —KAX; (x —x,) 





Utilizing the vectors from g(x), (4) and (5), a matrix is constructed. By 


investigating the determinant of this matrix, the controllability may be characterized. If 
the matrix is non-singular, then the system is weakly controllable. Conversely, if the 
determinant is singular (i.e. equals zero), then the characterization is undetermined and 


the Lie bracket method offers no insight to the controllability of the system. 





Compute the determinant: det| 8,[ 7.2], f[F.2]]]. 


M+x, eX, (x —x,)+ (A +«(x, + x,)) —KAx, (x,-x,)+ WA(A+«(x, +x,)) 
=> det|-(4+x,) «x,(x,-x,)—u(A+K«(x,+x,)) —KAx,(x,-2,)-HA(A+«(x,4+x,)) 
0 Ki eH; | —KAX,(x,—X,) 


pene (x op es =x )—M(A+K( x, +x;,))) | 


“(ats -x)-an (Ata +8) 


—A (cx, (x,-x,) + (A+ «(x +2,))) | 


—(-KAx, (x, —X;,) +pA(A+K(x, +x,))) 


soonest} 


= (ut x,)(«x, ( —x,))[2ua(A+«(x +x,))| 
+ (p+ x, (ex, (4 -% ))[-2uA (A+ K(x, +x,))| 


= Kx; (%,—%,)(%, —x,)[ 2wA(A+K(x, +x,))| 


=> det] a, 7.8]. FLA2]]| = 2AKx, (A+K(x, +x,))(x, -x,) 


A 
The determinant does not equal zero where x, 40, x,+x,=—-— and x,#x,. It 
K 


can be concluded that the system of the Phan-Thien-Tanner model under extensional flow 


satisfies Definition 3 and thus by Theorem 1, it is weakly controllable when x, #0 (or 
A A 

T,, #0) and x,+x,=-— (or 7,,+7,, #-—) and x, #x, (or 7, #7,,). This result can 
K K 

be summarized geometrically: 


Given R* = {(x,,x,,x,)|,x,,%, © R}, there exists three defined surfaces: 


S, = {(.%)5%5) 1% =0} 
S, = {(x,,%),%,)|4+«(x,+x,)=0} 
S 


3 meee aes =o} 


such that the system is weakly controllable at all points in R* \{S, US,US at : 
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B. HOMOGENOUS SHEAR FLOW 


A fluid in a homogeneous shear flow with rate y(t) is defined by velocity 
. : : : : 0 1 
v=(7(t) y,0), where the velocity gradient tensor is Vv = 7(t) a Gil Recall, a 


schematic plot of shear flow is given in Figure 2 of section I. This velocity is applied to 


the Phan-Thien-Tanner system (2), whereby the system’s components are investigated. 


=f -(Vv)T-T (Vv) +A +«(T)T =u] Vv +(Vv)" | 


=F-+0([) tI, saree n= artol[t Aue ‘| 
“lt POS Thee dpeoeemle ghee 


(ome 8 T, T. 0-1) (2h £ 
Siege oe =-(2+4(, 4) : a ul }4 a y(t) (6) 
T; Ty, T;, Thy 1 0 T, 0 














The system’s components are as follows: 


ST, =-(A+K(T,4+T))T + 277 (t) 
=>T,, =-(A+«(T, +Ty))T (7) 


For mathematical convenience, the following notion is introduced: 
qT, 
5 eae aa av om 
qT, 
System (6) can be rewritten as: 


dx 


nt (eRe 
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o 
a 
Ss 
—~ 
at 
— 
st 
| | 
——~ 
wm &d& 
+ + 
A A 
— 
coe! 
Ee ace 
Raleeta! 
—_ — 
_——_ > 
a 3 
———— es | 
Oy 
per 
at 
— 
Il 
rs | 
N 
j=) & 
———E—EE | 
= 
II 
R 
a, 
~ 
— 


i oe —(A+«(x,+2x,))x, 2, 
|r -(A+«(x,+x,))x, | 0 be 
qT, —(A+«(x,+x,))x, H+Xx, 


The following Lie Brackets are computed: Ez 3 | ; lf. Ez g ]| ‘ 


First Lie bracket: lf. 3 | =V 
2x, —(A+«(x,+x,))x, —(A+K«(x,+x,))x 2x. 
(78) 0 —(A+«(x,+x,))x, -~V —(A+«(x,+x,))x, 0 | 
—(A+K«(x,+%))% —(A+«(x,+x,))%, M+ X, 
0 2 —(A+«(x,+x,))x, 
ad x —(A+«(x,+x,))x, 
0 OM ie (A+K(x, +X,)) x; 
) 


0 

0 

1 
(A+x( (2x, +x, ) KX, 0 ON 
(A+x«(x,+2x,)) 0 | 0 


KX; (A+«(x,+%)) M+X, 


-2(A+«(x,+x,)) x, 2x, (4+«(2x,+%,)) 
ad [7.8 = 0 ap 2KX,X, 
( 2KX; 


- A+K(x,+x))x, 3 +(utx,)(A+K(x,+x,)) 


g-f-Vf-3 


0 
=| f.3| = 2KXyX; (8) 
2x; + (A+ K(x, + ;)) 


Second Lie bracket: EAERal =V| f.8|-f-Vi | #8 
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2KX,X3 





: —(A+K(x, +X,))X, 


; 0 0 | -(A+K(x +X,)) x, 


—(A+K(x, +x))% | 





(A+x«(2x,+x,)) KX, 0 0 
+ KX, (A +K(x, + 2x,)) 0 ; 2KX,X3 
KX, eo (A+«(x,+)) | 2Kx; + u(A+ K(x, +x,)) 


0 
-[F{F] eters} 25a | 


MX, + UX, +425 


2K xX, XpX; 
+ 2KX,x,(A+K(x, + 2x,)) 


2K xx; +(A4+K(x, +x, ))(2xx5 + (A+ K(x, + x,))) 
DEON 
=F 78]|= 0 (9) 
—2Kx;(A+Kx,)+ uA 
Construct a matrix with the vectors from g (x) , (8) and (9) to compute the determinant: 
det| #78], F[ 72]. 


2 
2h 0 DR Spa 





=> det} 0 2KX,X; 0 
M+x, 2x; +u(A+K(x,+x,)) —2Kx;(A+Kx,)+ pA 


=4Kx,x; | -2AKx,x, + PA-K°X,X, (u + x, ) - 2K x, | 





The determinant does not equal zero where x,4#0, x,#0, and 
2 MA 
KX,X)[uU+x,]+2x3|[A+Kx,]=—. It can be concluded that the system of the 
K 


Phan-Thien-Tanner model under shear flow satisfies Definition 3 and thus by Theorem 1, 


it is weakly controllable when x,#0 (or 7,,#0) and x,#0 (or 7,40) and 
A 
Kx a, [tay] +22 [At+nx]=4 (or KT, Typ [+Ty ]+ 272 [A+KT, ]=-—). This 
K K 
result can be summarized geometrically: 


Given R® = {(x,,x,,x,)|%,2,,%, © R}, there exists three defined surfaces: 


Si = (eee =0} 
S,= He sag a) ly =0} 
S 


Me {(4,.%2>%) | HA K(Kx,X, [w+x,|+2x5 [A +«x,]) =0} 
such that the system is weakly controllable at all points in R° \{S, CIS, es 3} : 


This correlates with the work of Renardy [1]. By manipulating the system from 


(7), it was found that the set of reachable states is given by: 


sa{(r, (1) n.(e)atale))-m2(,)> | 


e€@° -k 


where C is a constant given by: 


If 7,,(0)=0, then the reachable set includes the boundary point 
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The result obtained by Renardy is stronger than the one computed with the Lie 
Brackets in the sense that his result is on the global control whereas the later is on the 
local control. As mentioned earlier, weak controllability is practically more useful in 


reducing uncertainties and noises. 


C; SUMMARY OF THE PHAN-THIEN-TANNER RESULTS 


The Phan-Thien-Tanner model under both homogeneous extensional flow and 
shear flow, was shown to be weakly controllable on certain submanifolds via the Lie 
bracket method. Though Renardy [1] had already shown that the Phan-Thien-Tanner 
model under shear flow could be characterized, this thesis extended the concept by 
utilizing a different form of analysis, the Lie bracket method. Furthermore, the concept 
was broadened through the introduction and investigation of an additional flow rate, the 


homogeneous extensional flow. 


be 
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Il. JOHNSON-SEGALMAN MODEL 


The Johnson-Segalman model characterizes the behavior of non-Newtonian 
fluids, including special cases of Newtonian and Maxwell fluids. Additionally, it is a 


viscoelastic fluid model which was developed to allow nonaffline deformations [6]. 


This model has been used successfully and can be used to explain “spurt” 
phenomenon (a large increase in volume for a small increase in the driving pressure 


gradient, at a critical pressure gradient) [7]. 


The Johnson-Segalman model gives a viscoelastic constitutive equation and can 


be described by the following state equation: 


T-<—((vv)r +7 (vv)" )- (vey T+T(Vv))+AT = 2D (10) 
where the definition of its element are as follows: 


qT, Tp 


T. 


T : Stress tensor >T -| 
21 22 


(where 7,, =T7,, due to symmetry) 


T : Derivative of the stress tensor with respect to time 


a: Parameter describing polymer slip (a measure of the nonaffinity of the 
polymer deformation),where —1<a<1. Note: for the case where a =1, the 
model reduces to the Oldroyd-B model [8]. 


v: Velocity where v =(¥,,v,) 


A: Relaxation rate 


i: Elastic modulus 


Vy +(Vv)" 


D: Rate of deformation tensor (the symmetric part of Vv) > D= 5 


y: Control input (closely related to velocity) 
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A. 


HOMOGENOUS EXTENSIONAL FLOW 


A fluid in a homogeneous extensional flow with rate y(t) is defined by velocity 


y= y(t aes —y(t ) 2), where the velocity gradient tensor is Vv = 


2 


nak : 


This velocity is applied to the system (10), whereby the system’s components are 


investigated: 


=| a a HO es 
qT, 


fa Ts 


cal ale Tn 
lb ak: all 


N 


i 








_| Fi | -2O| Ti mh 
[Ti Ty) 2 T, —T,, qT, 

=| 4 Plana fu P hextolalg 
Tis [os 12 Ty 0 

_,| ti lan i ela 
| Tio 13 LT. T,, 0 





The system’s components are as follows: 
=> 1, =—AT, +(u+aT,,)7(¢) 


= 1, =-AT,, = (i+ aT) ¥(t) 


=>T,, + AT, =0 


nial 4) 
oer 


= u(Vv +(vv)') 


ae 
he eet 


Dee ch 


Fe )fteenetemnean 


= 2 eae | 
een gol, {i 


(1) 


(12) 
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Note that Equation (12) can be solved exactly: 
=T,, =T,,(0)e" 


Consequently, the Johnson-Segalman model is not controllable. Nonetheless, the state 


space has a stable invariant subspace T,, =0. It is this subspace that all trajectories of the 
system, under any control input, asymptotically move toward the subspace of 7,,=0. As 


a result, the decisive behavior of the control system is characterized by the reduced 


system on the stable subspace, 


(13) 


For mathematical convenience, the following notion is introduced: 


tbe 


Then system (13) can be rewritten as: 


ole -( || M+al,, \ro 
T,, | [-ATy | |-(“+aTy) 
The following Lie Bracket is computed: Ee 3 | : 


=| f.e]-VE-f-VF-3 


Oa 
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= 1 Pl yee Mira -| eal ais 


A(u+ ax, ) adx, —A(p+ax,) —Au 


Construct a matrix with the vectors from g(x) and (14) to compute the determinant: 
se(e[7.a]] 

= A 
soa[a[7a]]=oe fatale) 


The determinant does not equal zero where x,#x,. It can be concluded that 


since the sub system (13) of the Johnson-Segalman model under extensional flow 


satisfies Definition 3 and thus by Theorem 1, it is weakly controllable when x, # x, (or 


T,, #T,, ). This result can be summarized geometrically: 
Given R® = 1a) lx,,x,€ Ri , there exists one defined surface: 
S,= 1(2.35)) Ix= x, 


such that the system is weakly controllable at all points in R? \{s i} : 


B. HOMOGENOUS SHEAR FLOW 


A fluid in a homogeneous shear flow with rate 7(t) is defined by velocity 
; spe er 
v= ( y(t) y,0) , where the velocity gradient tensor is Vv = y(t) Ae This velocity is 


applied to the system (10), whereby the system’s components are investigated: 
qT T,, 2 0 0 qT Ty T, Ty 1 0 
(rel) le ack F baole eae | 
2 1 0 T T,, qT) Ty 0 0 qT T,, 


= (Vv +(vv)') 
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1 
u+5((4-I)T+(4+)T) (a-1)T, 
The system’s components are as follows: 


=> T,, =-AT,,+(a+1)T,,7(t) 
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=> T,, = -AT,, +(a-1)T,,7(t) (16) 


=f, =-AT, +] een, +(a+IT) [70 


For mathematical convenience, the following notion is introduced: 


System (15) can be rewritten as: 


<= F (3) +8(8)u 
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—AXx, (a+1)x, 
Where jets} g(x) (a-1)x, 
—Ax; 


qT, —AX, 
=>|T, |=| -Ax, |+ 
qT, —AX; 








(a+1)x, 


w+ ((a-1) +(a +1)x,) 


w+5((a-1) x, +(a+1) 5) 








v(t) 





We will now compute the following Lie Brackets: Re g | : EARZ g ] j 


Compute first Lie bracket: lf. 3 =Vg: jf —Vf-g. 


=| f.]=V 








(a+ 


1) —-AX, 
(a-1)x, fe 


w+s((a-t) +(a+1)x,) 








(a+1)x, 
(a —1)x, 


w+ 5((a-1) x +(a+1)x,) 


0 0 (a+1) | | -Ax, 
0 0 (a-1) fe 
ae), oe 
-A 0 0O | (a+l)x, 
| —A | (a-1)x, 
0 O -A 
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w+ 5 (at) +(a+1)x,) 





—x,(a+1) x,(a+1) 
=| f.g|=4 —x;(a-1) +A x;(a-1) 
-(4(a-1)+2(a+1)] w+5((a-1)x,+(a+1)a) 
0 
=[F.8]=| 0 oy 
LA 
Compute second Lie bracket: |F[f2]| =V| f.8|-f -Vf+[ f.8 | 
0 | | -Ax, =AX; 0 0 
=(A[f.8]]=V] 0 |-|-4% |-v)-a~ |] 0 | =[F[Aa]]=| 0 (18) 
HA | | -Ax, —AX, | | HA wa? 


Construct a matrix with the vectors from g(x) , (17) and (18) to compute the 





determinant: det| 8, 7.8], F[ 72]. 


(a+1)x, 0 
=> det (a-1)x, 0 0 |=0 


w+5((a-1)x +(a+1)x,) UA pA? 








The determinant equals zero everywhere and thus is singular. It can be concluded 


that the controllability of the Johnson-Segalman model under shear flow characterization 


cannot be determined by g, Ear and Besa Nevertheless, it is possible to 


obtain control invariant manifolds in the state space that are weakly controllable. This 


needs more advanced geometric analysis and it is beyond the scope of this thesis. 


It should be pointed out that by some clever direct analysis, Renardy [1] was able 


to derive some complicated expressions for the reachable set. His result can be 


ZS 


summarized as follows. We first introduce a new variable which is a linear combination 


of the two normal differences T,, and T,,: 
a 1 
Z = 5 (hi + Tn )+5 (En =i) : 


12° 


Then the Johnson-Segalman system can be reduced into a system of ODEs for (7,,,Z): 


Ti, =-AT,, +(Z + #) y(t) 
Z =-AZ +(a’ -1)T,,7(t) 


Using arguments analogous to that for the Phan-Thien-Tanner model, Renardy derived 


the reachable set S defined in the (7,,,Z) plane: 


12? 
S ={(n. (1,).Z(t,)) 
Here in the inequalities the lower and upper bounds are solutions of the IVPs: 
d 
L(t) =-224 (0) +24 -1e+ 240)" | 
d 
= 219, (1) + An| -n—)20s (1)? 


so-so z0 920) eV) 


Ali, )s4Z(1,) +a2(1,)+5(-@ nal) <A(i)| 





If 7,,(0)=0, then the reachable set includes the boundary point 


(Tia(t,).2(t,)) =(0.-1+ 2a (t,)+ 4°) oF (Tia(t,).2(t,)) =(0.-1- 20h (0, #2"). 


Again, Renardy’s result is mathematically stronger than the one computed from the Lie 


bracket method whereas the later is more practical. 


Cc. SUMMARY OF THE JOHNSON-SEGALMAN RESULTS 


The Johnson-Segalman model with homogeneous extensional flow showed to be 
weakly controllable for a subsystem, whereas, the shear flow was undetermined. Though 


Renardy [1] had already shown that the Johnson-Segalman model under shear flow could 
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be characterized, this thesis extended the concept by utilizing a different form of analysis, 
the Lie bracket method. Furthermore, the concept was broadened through the 
introduction and investigation of an additional flow rate, the homogeneous extensional 


flow. 


zi) 


THIS PAGE INTENTIONALLY LEFT BLANK 


28 


IV. DOI MODEL 


The Doi model for rod-like molecules is known for its ability to represent the 
characteristics of liquid crystal polymers in a solvent. The single molecule position- 
orientation distribution function is a fundamental element of the model. Exchanges 
between molecules are represented by a mean-field potential and the model can represent 
mean-field kinetic theory. Rods are subject to Brownian force due to the fact that they 
interact with other rods and with the flow. Generally, the model is a microscopic 
Fokker—Planck type equation (Smoluchowski equation depicting the convection, rotation 
and diffusion of the rods) coupled with a macroscopic Stokes equation (pertaining to the 


hydrodynamics) [9, 10]. 


A fundamental attribute of the Doi model is that it has the capability to describe 
both the isotropic and nematic phases. If the interaction strength of the rods is effectively 
strong or if the rod concentration is high, then the system tends towards a nematic phase 
and the rods have a propensity to line up. If the rods randomly orientated, then the 


system is an isotropic phase [9, 10]. 


The full Doi orientation tensor theory is developed after the microscopic Fokker- 
Planck equation is projected onto a second-moment description using closure rules. The 
major ingredient in this tensor theory is the second-moment tensor which describes the 
orientational distribution of the ensemble of rod-like macromolecules. The orientation 
tensor is traceless and symmetric. It is the basis for micron-scale light scattering 
measurements of primary axes (directors), degrees of molecular alignment 


(birefringence), and normal and shear stress measurements. 


This section is devoted to the controllability of two-dimensional Doi tensor 
model. The study of two-dimensional liquid is physically motivated by monolayer films. 
As a result of their stability and nonlinear optical characteristics, thin films of liquid 
crystalline polymers have drawn a significant amount of attention contrast to materials of 
low molecular weight. The controllability of this model under homogeneous flows will 


be explored for the first time in the literature. Since the bulk properties of liquid 
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crystalline polymers (LCPs) depend on the microstructure, especially the orientation of 
the molecules, the ability to control orientation tensor will open the door for achieving 


the desired strength and materials properties of LCPs. 


The two-dimensional Doi model for liquid crystalline polymer is described by the 


following state equation: 
0=20-09+ a[D@Q +QD] +aD—2aD:Q [o Z| —~6D°F (Q) (19) 


where the definition of its element are as follows: 


2, 1 0, 


0, a (where —Q,, =Q,, and Q,, =Q,,) 


Q: Orientation tensor , > Q = 


Q : Derivative of the orientation tensor with respect to time 
oe 1 z 
Q:: Vorticity tensor > QO = ae —Vv ] 


v: Velocity where v =(v,,v,) 


a: Dimensionless parameter depending on the molecular aspect ratio, 


D: Rate-of-strain tensor => D = sly + vy" | 


F(@)-(1-XJo- no" +no:0[ 0+) 


QO, -Q,, 0 (0,,) +(Q,) 
+2N ((Q,.) +(Q,,) 1 
Q,, Lat 


N: Dimensionless concentration of nematic polymers, 


D°: Averaged rotary diffusivity or relaxation rate 
30 


Notation [oO :Q| = tr(QQ" ) 


A. HOMOGENOUS EXTENSIONAL FLOW 


A fluid in a homogeneous extensional flow with rate y(t) is defined by velocity 


/ 1 O 
v= [70 : ,-7(t) | , where the velocity gradient tensor is Vv = mek i z 


This velocity is applied to the system (19), whereby the system’s components are 
investigated: 


= 9-020)! ale I A Al, | 7(t) k | 
2 Oe QO, -O,, Q,, Q, Oi] 
. 1 Q, QO Q, A 11 0 _ En? 
ely lle: SiMe: “Silo i] re 
= 9-20) & oe || Pale zal i 
2 —O,, Q,, 0, QO, 2 es 
. Q, A Q, QA id 0 0 
onl le oF 20 || 90°10) 
1 0 . Oy, Q, Eb 
va * Fer e(2a.)| 0, va | 


coh eels 0 | 
2 Q,, =O 0 : 





Q,+Q,, 
_6)N° 1 
OD, Q, ss 0, 
2 2 
+2N(Q;,+Q;) 
Q,, -Q,, a 
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[1-2/6 eae 0 | 
2 2, -O, 0 Q;, +0), 


=O =0D, On 0; 
+2N(Q2+02)|  * 
2, Lit 
0, 0 a 1 0 0, 0, 1 1 0 . 
ick a al eae) Ae re 
[1-2/2 Q, -w{e +Q7, 0 | 
2 2, -O, 0 Q,,+ Or, 
aoe Qt O% 
+2N(Q? +03) 2 
Q,, =i (20) 
2 
ie 0 ee: G, 
ral] *  f-29,) * Wo 
0 Qu-5 0, eae 


We get the following two component equations: 


= O,, = —6D? ([-T}o, +2NQ,, (Q7, +03) jra{ 3-203 }(0) 


: N : 
> 0, =-6D"O, ((-%) 2N(O? +03) |-240,0:7(0) 
Regard y(t) as a control input and consider the problem (19) with initial state Q, 
and final state Q,. The system is controllable if for any choices of Q, and Q, there exists 


a control y(t) on the interval (0,7). 


Using the nematic relaxation time scale Dp’ the flow field and orientation 


r 


dynamics of (19) can be non-dimensionalized. The key dimensionless parameters are 
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then the Peclet number Pe(t) = I (the shear rate normalized with respect to nematic 


r 


relaxation rate) and the dimensionless concentration parameter N . 


= O,, =D? -6g,,{1-+2N(0%, +03) }+a{ 5-20; ) P(r) 


o 





= Q,, =D, | -6Q,, a-24 2N (Qi, £05 )}-240,0,0e(0) 


Applying t =1D°, the equations become dimensionless and the nematodynamic 


model (19) can be expressed as: 


= 0, =-69,(1-+20 (9% +03) +a{ 2-208 ]Pe(0 


2 
< (21) 
=> 0,, =-60,, [ arr: 2N (Q7, +Q7, ) —2aQ,,0,,Pe(t) 
For mathematical convenience, the following notion is introduced: 
x= = 
xy 2, 
Then system (21) can be written as: 
OE. Fp 
x u 
B= F(8)+8(3) 
-65,[1-24 +2N( (x; +X5 ) [ " 
2 eke Wal ek 
Where (x)= - g(x)=| \2 , u=Pe(t) 
—6x, [i- a 2N ( (x; + 5) Nf —2ax,X, 


Compute the following Lie Bracket: Fe 3 =Ve-f—Vf-g. 
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7 [ar 0 -65,(1-L4-2n (x? +92) 
[F.8]= —2ax, eel 


N 
1-2 2N (3x7 +x) —24Nx,x, of 5-237] 


—24Nx,x, —6 [1 = = +2N (x; + 3x5 ) —2ax,X, 


—6x, [1S 2m (ai +32) 


24a (1-4 20 (x5 +x) 
Fa-| 
24ax,x, -34 2N (x; + «) 

—6a [ -$+ 2N (3x7 +X5 iG - 2x + 48aNx, x; 


1 
—24aNx,x, [5— 2a J+ 120K [34 2N (x; +332) 


N N 1 

6a; (1-42 (x) +49) }+(1- +20 (aa? +22) [5-297 )-snay| 
N I N 

12a 2(1- +20 (ap +2) +2n(5 237] [ v+20(x3+32))| 


6a ang 1-L ann? ]a{1-L 4-20 (397 +28) |[ 3-209 
2 2 2 


-2N (x; +x3)} 





| a | 
a) 
oa 
Ld 
ll 
w|= 


12ax x, f + 


6a{2x{(1-L-an (x? +23) )+3-“aw(aaf +2) 
| FB |= (22) 
12ax,x, {+2 an (af +x8)} 





Construct a matrix with the vectors from g(x) and (22) to compute the determinant: 


se( [7-4] 
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1 
a er 6a4 2x, 1-2 (af 4x8) joan (39 +) 
2 2 2 4 
=> det 
—2ax,X, 12ax,x, {i+ -an (af +22) 


=124° xx, [5-27 f+ Lan (9 +x5 ) 


N 1 oN 
vi2a?ss {28 (1--an (+a) +2 (ae) 
1f{. N 1 oN 
=12a’x,x, {i+ San (at +22) }-ang +d Ma w(aat +d) 


> det| 8] f, z|| =12a?x,x, 


The determinant equals zero where x,#0 and x, #0. It can be concluded that 


the controllability of the sub system from the Doi model under extensional flow satisfies 


Definition 3 and thus by Theorem 1, it is weakly controllable when x, #0 (or Q,, #0) 


and x, #0 (or Q,, #0). This result can be summarized geometrically: 


Given R* = {(%4.%) Ix,,x, € R} , there exists two defined surfaces: 


such that the system is weakly controllable at all points in R* \{s ,VUS ot ; 


B. HOMOGENOUS SHEAR FLOW 


A fluid in a homogeneous shear flow with rate y(t) is defined by velocity 


; 1 
v= (7(t) y,0), where the velocity gradient tensor is Vv = (09 i 
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This velocity is applied to the system (19), whereby the system’s components are 


investigated: 
29-20" 2: 2. |_| alk il 
“1 OjLQ. Or Oe SOe al +0 
at ile Q,, He 05 lf Nita O ‘| 
2 1 0 0. -O, 0, =075 1 O 2 [1 0 
y a 2, Q,, Q,, Gi 1 1 0 | 
—ay(t) 1 0 : 0 0 5% ts a 


29-10) QO, —O,, -| > Q, al? —O, +| & Q, 





-~6D°F (Q) 


— 






































2 |L-O, -Q.} LQ, @ 1) 2 (LQ, Qs} LO, |, 
OT ob ereve(e Si feale i) 
=9-H0) “oo fea Ce : aX! i 
2 |-20, -20,, | DIO) 20s 2 |1 0] 
-247(1)04{| 0" Pace i] erro) 
































0, “O's 2. Yok~ Oy 
Q, OA 1j}1 0 
-2«7(N0,{|6 O A 1 
(#2 On ]_y[(2ir@) 9 | 
2)|Q. Qi 0 (Q;, +07) | 
—6D° 1 
Qt Q, 

+2N(Q;, +7.) 
QO, cae 
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[1-2/2 2, La 
2/19. —Qy] 0 








(OF +@5 ) 0 











(O, +90), 
)=-6D° 1 
ne ’ Qt Q,, 
+2N(Q;, +07) 
0, ts. 
Q,, -O,, | ks 0 a F 1 
+a +— 
-O, -Q,, 0 Q,, 2{1 0 ; 
Bact Geen) 
ie Q,, | ) 
2aQ,, +— 
Q,, 2), | 2|0 1] 


The system’s components are as follows: 








O,, =-60'0,[1-+2N(Q} +03) |+0 (1-2aQ,,)7(t) 


N 


Or, = -6D'0,, a-24 2N (Or +Q )}+{-2, ‘5 200% (0) 


(23) 


Regard 7(t) as a control input and consider the problem (19) with initial state Q, 


and final state Q,. The system is controllable if for any choices of Q, and Q, there exists 


a control y(t) on the interval (0,7). 


Using the nematic relaxation time scale 


—, the flow field and orientation 


r 


dynamics of (19) can be non-dimensionalized. The key dimensionless parameters are 


then the Peclet number Pe(t) = ot (the shear rate normalized with respect to nematic 


r 


relaxation rate) and the dimensionless concentration parameter N . 


N 


=> Q,, =D? +60, f-34 2N((Q,) +(02))}+ Q,, (1-2aQ,,) Pe(t) 
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: N 2 2 a 2 | 
> 2, = D: +60, a-34 2N((Q,:) +(Q.) )}+{-2, +5 ~2a(Qn) Jre(o 
Applying t =tD°, the equations become dimensionless and the nematodynamic 


model (23) can be expressed as: 


=> 0, =-60,, a-34 2N((Q,) +(0.))}} Q,, (1-2aQ,,) Pe(t) 


(24) 
. N 
> Q,, = —62,, a-24 2N((O,) +(Q2)))}*(-2: o 2a(Q..) Pe) 
For mathematical convenience, the following notion is introduced: 
{He 
x= = 
Xy Q,, 
Then system (24) can be written as 
& =F (3)+8(8)u 
N 2, 2))| 
—6x, - Zs 2m(s +X; ) x, (1-2ax, ) 
Where f(X)= , a(x)= pea u= Pe(t) 
-61,{1-42(x¢ +32)] i ae 


Q, 





-{2]- -645[ 1-24-26 (sf+32) | x, (1-2ax,) Jt 


-61,(1-L 420 (x3 +32) 


II 


Compute the following Lie Bracket where lf. 3 | Vg: f Vp g. 
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—6x, a-24 2N (x; +32) 
—6x, [1-2 2N (9 +32) 


-6(1 3+ 2N (3x7 +x5 ) —24Nx,x, | x, (1—2ax, ) | 


a 2 
~24Nx,x, -6[1-+2N (x8 +3%)) Mr Se 


N 
—6x, f-34 2N (x +32) ][l-4a5] 


= [Fa i N 2 2, 2 
6(1-S 420 (x xg) [n+ dan] 
—6x, f-34 2N (3x; +X5 )Ja- 2ax,)+ 24.Nx,x, [s at Dax: 
—24Nx,x5 (1-205) +6{1-24 2N (x; +3x5 )[ =e das’ 
12ax,x, [i+ aw (a? +2) 
=> F:8|= (25) 


sa(4u 14a | 2-203 (9 ta)+o(¥ +2)) 
Construct a matrix with the vectors from g(x) and (25) to compute the determinant: 
det| 3, 7.8]. 
1-2 12 144 on (242 
ai ax, ) aX,X, | 1+ 5 (x; +x;) 
=> det i ; 
a 
—%x, +5 2ax, sa{4u ei4an| 220 (9 +a5)+5(a +22)) 
= x,(1-2ax,)3af 422-4144] == 202 (2 +22) +(22 +22) 
8 2 
—12ax,x, 1+ —2n (af +22) [[-1+$-2ax8] 
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(1- 2ax, ) (4x3 sh das) -5~ 243 (x9 +38) +5 (97 +22) 
= 3ax, . 
- 4x( a 5 —2ax; }: | -2.Nx, +8Nx, (x; +X5 )| [-s os ae Dax 





(1-2ax,)(4x? +1) an( x45 Zax: 


= 3ax, -2-4e (x; +x3)+(a + 2x; )-2ax, [-t-4(x xg) +(af-+2a2) 


EX |-1+ A(x; +x5 V[-« i: 2a; 


=3ax, Jee (a +25 )—4x7 (a7 +25) + (27 +37 +225) 
-tag[ 1 auc(af ed) a(ul 2e) a(t) eae (+2) 


x +37) )+1- 4ax, 


129] [-t-4[d eae +8) 2(0 +8) 


= 3ax, | 5-8 x; 2+ x2) +4(sf 438) |e4(af taf) 1d 


| 
; 
| 
al 


=>det f.8 ||= Sa, |4( 2+ 92)(-2M (9-42) 4 41)-N ada 


The determinant equals Zero where x, #0 and 
2 2 2 2 N 
A(x, ay )(-20v (x; + x5 ) +N+ 1) ae +1—4ax,#0. It can be concluded that the 


controllability of the sub system from the Doi model under shear flow satisfies 


Definition 3 and thus by Theorem 1, it is weakly controllable when x, #0 (or Q,, 0) 


and 4(.x; +2)(-2N (x; bad) + N41)-2 41-day #0 


40 


(or 4(Q;, + Q}, )(-2N (Q;, +Q;, ) +N+ 1) - > +1—4aQ,, #0). This result can be 
summarized geometrically: 


Given R* = 1( 4535) Ix,,x, € R} , there exists two defined surfaces: 


S, 1 (ayaa) =0} 


5, ={(5.a)l4(x+28)(-2N(ah +x) rv 41) stax =o] 
such that the system is weakly controllable at all points in R \{s ae) S,} : 


C. SUMMARY OF THE DOI MODEL RESULTS 


The sub system from the Doi model under both the homogeneous extensional 


flow and shear flow showed to be weakly controllable via the Lie bracket method. 


4] 
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V. SUMMARY OF MODELS IN TABLE 


The following table depicts a snapshot of the results found in this thesis. 

















MODEL Homogeneous Homogeneous 
Extensional Flow Shear Flow 
Phan- weakly controllable weakly controllable when T7,, #0 and 7,, #0 
Thien- when 7,,#+7,, and ‘i A 
Tanner Tis anh 1 20 and wT Tn [M+T]+ 20 [4 + xT,,] =" — 
Johnson- weakly controllable undetermined with the chosen vectors 
Segalman | when T,, #7,, 
Doi weakly controllable weakly controllable when @Q,#0 and 
when Q,, #0 and N-2 


Q,, #0 








(Q/, +O; )(-2N (Of, +O3)+N +1)-a0,, 4 





Table 1. Constitutive Model Overview 
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VI. CONCLUDING REMARKS AND FUTURE WORK 


The work covered in this thesis investigated the controllability of three popular 
constitutive models under homogeneous extensional and shear flows via the Lie bracket 
method. By applying tools from the geometric control theory, submanifolds on which 
each system or its subsystem is weakly controllable were explicitly identified. 
Characterizing a system by means of this method is a fairly new concept and it has 


applicability to a wider variety of models. 


Potential extensions of this paper include the investigation of controllability in 
three-dimensional systems, though it definitely becomes quite challenging. Note that this 
thesis only gave weakly controllable set. Identifying reachable set for each model under 
extensional flow should be a future work. Another possibility for future work is the 
exploration of controllability of the state equation under nonhomogeneous flows. This 


requires analysis of PDE systems and it is undoubtedly more complicated. 


45 


THIS PAGE INTENTIONALLY LEFT BLANK 


46 


LIST OF REFERENCES 


[1] M. Renardy, Shear Flow of Viscoelastic Fluids as a Control Problem, Journal 
Non-Newtonian Fluid Mechanics. 131, 2005, pp. 59-63. 


[2] M. Renardy, On Control of Shear Flow of an Upper Convected Maxwell Fluid, 
Z. Angew, Math. Mech. 87, 2007, pp. 213-218. 


[3] A. Isidori, Nonlinear Control Systems, 3rd edition, Springer, 1995. 


[4] H. Zhou, W. Kang, A. Krener, Steady Flow Field Effects On The Control Of 
Maxwell Materials, Preprint 


[5] H. Giesekus, A Simple Constitutive Equation for Polymer Fluids Based 
on the Concept of Deformation-Dependent Tensorial Mobility, Journal 
Non-Newtonian Fluid Mechanics. 11, 1982, pp. 69-109. 


[6] M. Johnson Jr. and D. Segalman, A Model for Viscoelastic Fluid Behavior 
which allows nonaffine deformation, Journal Non-Newtonian Fluid Mech. 
3, 1977, pp. 255-270. 


[7] M. Elshahed and M. Haroun, Peristaltic Transport of Johnson-Segalman Fluid 
Under Effect of a Magnetic Field, Mathematical Problems in Engineering 6, 2005, 
pp. 663-677. 


[8] V. Ervin, J. Howell, and H. Lee, A Two-Parameter Defect-Correction Method for 
Computation of Steady-State Viscoelastic Fluid Flow, AMS Mathematics subject 
classifications, 65N30, pp. 3. 


[9] M. Doi and S. F. Edwards, The Theory of Polymer Dynamics, Oxford University 
Press, 1986. 


[10] CC. Luo, H. Zhang, and P. Zhang, The Structure of Equilibrium Solutions of the 
One-Dimensional Doi Equation, LMAM and School of Mathematical Sciences, 
Peking University, Beijing 100871, People’s Republic of China, 2004, pp. 379. 


[11] R.B. Bird, R. C. Armstrong and O. Hassager, Dynamics of Polymeric Liquids, on 
edition, Wiley, 1987. 


[12] R.G. Larson, The Structure and Rheology of Complex Fluids, Oxford, 1998. 
[13] H.Nijmeyer and A. J. van der Schaft, Nonlinear Dynamical Control Systems, 


Springer, 1990. 


47 


THIS PAGE INTENTIONALLY LEFT BLANK 


48 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
Ft. Belvoir, Virginia 


Dudley Knox Library 


Naval Postgraduate School 
Monterey, California 


49 


